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Abstract 

Based on an approach recently proposed by us, simultaneous x 2 - anar y ses are per- 
formed for elastic scattering, direct reaction (DR) and fusion cross sections data for the 
6 He+ 209 Bi system at near-Coulomb-barrier energies to determine the parameters of the 
polarization potential consisting of DR and fusion parts. We show that the data are well 
reproduced by the resultant potential, which also satisfies the proper dispersion relation. 
A discussion is given of the nature of the threshold anomaly seen in the potential. 
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A great deal of effort has recently been focused on studies of the so-called threshold 
anomaly pi]2] (rapid energy variation in the strength of the optical potential) in heavy 
ion scattering induced, particularly, by very loosely bound projectiles such as 6 He [[|, 
6 Li and 9 Be The experimental results accumulated so far indicate that the 

imaginary part of the optical potential, W(r;E), extracted by analysis of the elastic 
scattering data, does not show such an anomaly as is observed in the potentials for 
normal, tightly bound projectiles. For tightly bound projectiles, W(r; E) at around the 
strong absorption radius r = R sa is found to decrease rapidly as the incident energy E 
falls below the Coulomb-barrier energy E c , and eventually vanishes at some threshold 
energy E . Contrary to this, for loosely bound projectiles W(R sa ; E) remains large at 
energies even below E c |||5||||^] . 

The reason for W(R sa ; E) being so large at low energies has been ascribed to the 
weak binding of the extra neutrons to the core nucleus, leading to breakup. In fact, 
the breakup cross sections have been measured for these projectiles P- pTjl , confirming 
that they are indeed large, even larger than the fusion cross sections at E ~ E c . It 
was argued that since the energy dependence of the polarization potential due to the 
breakup must be weak, one might not be able to observe a noticeable energy variation 
in the resultant potential when the breakup cross section is larger than the fusion cross 
section as for loosely bound projectiles. 

The threshold anomaly of W(r; E) observed for tightly bound projectiles may be 
ascribed to the coupling of the elastic and fusion channels |TJ|]- This is substantiated by 
the fact that the threshold energy E of W(r;E) (i.e., the energy where W(r,E ) = 0) 
agrees very well with that of the fusion cross section up, or more precisely the threshold 



energy of S(E) = \JEop [|13|]. It is thus natural that if the breakup cross section is larger 
than the fusion cross section, and if one is concerned only with the total W(r;E), the 
rapid change in the fusion cross section and the anomaly would not show up clearly in 
the total W(r;E). 



2 



Insight into this problem may be obtained if one decomposes the total W (r; E) into 
the direct reaction (DR) and fusion parts, Wo(r;E) and W F (r;E), respectively, and 
determines them separately [H. The aim of the present study is to make such a deter- 
mination of Wd(t;E) and W F (r; E) by performing simultaneous x 2_ana ly ses of elastic 
scattering, DR (breakup), and fusion cross section data. We take the 6 He+ 209 Bi system, 
for which data are available not only for elastic scattering ||, but also for breakup || 
and for fusion |15| . Following Ref . || , we identify the breakup cross section with the DR 
cross section. Optical model analyses of the elastic scattering and total reaction cross 
section data have already been presented in Refs. |3|.[16| . The present analysis is thus an 
extension of the previous studies. 

The optical potential U we use has the following form; 

U = U c (r) - [V (r) + V(r; E) + iW(r; E)], (1) 

where Uo(r) is the Coulomb potential, whose radius parameter is fixed at a standard value 
of r c =1.25 fm, and Vo(r) is the energy independent Hartree-Fock part of the potential, 
while V(r; E) and W(r; E) are, respectively, real and imaginary parts of the so-called 
polarization potential |T7[| that originates from couplings to reaction channels. W(r\ E) is 
assumed to have a volume-type fusion and a surface-derivative-type DR part. Explicitly, 
Vo(r) and W(r;E) are given, respectively, by 

V (r) = V f(Xo) (2) 

and 

W(r; E) = W F (r; E) + W D (r; E) = W F (E)f(X F ) + AW D (E)a D ^§^-, (3) 

where f(Xi) = [1 + exp(Xj)] -1 , with X,- L = (r — Ri)/cti (i = 0, D and F), is the usual 
Woods-Saxon function. The real part of the polarization potential is also assumed to have 
DR and fusion parts; V(r; E) = V F (r; E) + Voir; E). Each real part may be generated 
from the corresponding imaginary potential by using the dispersion relation pf; 



E-E r°° W(r- E') 

V t (r; E) = K-(r; E s ) + ±^-^P / dE>- ' -, (4) 

K 3 J lV ' sj 7T Jo (E' — E s ) (E' — E) yi 

where P stands for the principal value and Vi(r; E s ) is the value of the potential at a 

reference energy E = E s . Later, we will use Eq. (4) to generate the final real polarization 

potentials Vf(t;E) and Vd(t;E), after Wf{t]E) and Wd(t;E) have been fixed from 

X 2 -analyses. For Vo(r), we simply use the potential determined for the a+ 209 Bi system 



at E=22 MeV [pLS[] , assuming that all the unusual features of the scattering may be 
described by the polarization part of the potential, particularly by the DR part. The 
parameters used for Vo(r) are Vo=100.4 MeV, ro=1.106 fm, and ao=0.54 fm. 

The unusual behavior of the elastic scattering and DR data for loosely bound projec- 
tiles can most dramatically be seen in plots of the ratios of the elastic differential cross 
section (daE/dfl), and the DR cross section (dap,/dQ), to the Rutherford scattering cross 
section (da c /dQ), i.e., 



as a function of the distance of closest approach D (or the reduced distance d) |T9| , pO 
that is related to the scattering angle 9 by 

D = i( AT + ^) = \o a (l + -^J with D = (6, 

Here D is the distance of closest approach in a head-on collision (s-wave). Further, 
(Ai, Zi) and (A 2 , Z 2 ) are the mass and charge of the projectile and target ions, respec- 
tively, and E is the incident energy in the center-of-mass system. 

In Fig. 1, we present such plots for two incident energies of £"=18.5 and 21.9 MeV 0. 
As seen, P E is close to unity for large d, but starts to decrease at an unusually large 
distance of d =2.2 fm (= dj, interaction distance). This value is much larger than the 
usual value of d\ «1.6 fm for normal, tightly bound projectiles. On the other hand, it 
is remarkable that the sum, Pg + Pjj, remains close to unity until d becomes as small 
as ~ 1.7 fm, implying that the absorption in the elastic channel up to this distance, and 
the unusual character of the scattering data, is due to the breakup. 



Since the theoretical cross sections are not very sensitive to the real polarization 
potential, we tentatively treat it in a rather crude way in carrying out x 2_ana ly ses j 
we simply assume Vi(r; E) has the same radial shape as the imaginary part Wi(r;E): 
Vi{r;E) = Vi(E)(Wi(r; E) /Wi(E)), Vi(E) being the strength of the real potential. We 
then carry out x 2 -analyses treating Wf{E) and rjj as adjustable parameters, keeping 
all other parameters fixed as Vp=3.0 MeV, r>=1.40 fm, ap=0.55 fm, Vd=0.25 MeV, 
14^=0.40 MeV and od=1.25 fm. The necessity of varying ao or rjj as a function of 
E has been shown in previous studies [p|,|l6|, and in the present work we take rjj as 
a variable parameter to study as a function of E. In the x 2 -analyses, data for elastic 
scattering, angle-integrated total DR, and fusion cross sections at E =14.3, 15.8, 17.3, 
18.6, and 21.4 MeV are employed. 

The values of Wp(E) and vd(E) fixed from the x 2_ana lyses are presented in Fig. 2 
by the open and the solid circles, respectively. Each set of circles can be well represented 
by (in MeV and fm, respectively, for Wf(E) and ro{E)) 



W F {E) 





1.25(^-15.4) 
4.0 



and 



r D {E) 



1.73 



E <15.4 
15.4< E <18.5 
18.5< E 

, E <14.0 



(7) 



1.73 - 0.03(^-14.0) , 14. 0< E <21.4 
1.508 , 21.4< E. 



(8) 



Note that the threshold energy E =15A MeV, at which Wp(E) = 0, is set equal to that 



of the linear representation of quantity S(E) = s/Eof oc (E — E ) discussed earlier. 



Kolata et al. |~5| found the value to be 15.4 MeV, which is used in Eq. (7). At this 
moment, we have no experimental information on r^-values below 14.0 MeV and above 
21.4 MeV. Thus, in Eq. (8), we tentatively set td(E) to be a constant as 1.73 fm for 
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E <14.0 fm and 1.508 fm for E >21.4 MeV. Note that the values of r D (E) at £=18.6 
and 21.4 MeV agree well with those determined by Mohr ||T6f . 

Eqs. (7) and (8), together with other parameters used for Wf(t;E) and Wd(t;E) 
as mentioned above, completely fix their values in the energy range between £=14.0 
and 21.4 MeV. In order to display the energy dependence of the potentials, we present 
in the lower panel of Fig. 3 the values of Wp{r; E), Wc(r; E), and the sum W(r; E) = 
Wf{v; E) + Wn(r; E) at a strong absorption radius r = R sa = 13.0 fm. It is remark- 
able that Wp{R S a] E) plotted by the dotted line exhibits a threshold anomaly (strong 
energy variation) similar to that observed for tightly bound projectiles. However, the 
Wo{Rsa] £)-values are rather flat as a function of E and have a magnitude much larger 
(by about a factor of 5) than the values of Wp{R sa \ E). Therefore, the threshold anomaly 
in Wf(R S cl', E) does not manifest itself in the total W(R sa ] E) plotted by the solid line. 

In order to generate the real part of the polarization potential by using dispersion 
relations, we need to know the imaginary potential in the entire range of E. Eq. (7) with 
Of=0.55 fm and rp=1.40 fm is enough for calculating Wf(t;E) in the entire £-range. 
For the fusion potential, since the geometrical parameters are energy-independent, the 
dispersion relation is reduced to that for the strength parameters Vf(E) and Wf(E), and 
the closed form for the expression has already been obtained |[| as 

V F (E) = V S (E S ) + -W F (E b )[e b ln\e b \ - ejn\e a \), (9) 
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where 

(E-E a ) (E — E b ) 

e a = j- and e b = j- (10) 

[E b — E a ) [E b — E a ) 

with E a =15A MeV and E b =18.5 MeV. The value V S (E S ) chosen is 3.0 MeV at £,=18.5 
MeV. 

For Wd(t;E), some care must be taken with the magnitude. To do the initial x 2 - 
analyses, use was made of Wd=0A MeV with od=1.25 fm in fixing the r^(£)-values 
given by Eq. (8). The constant value of Wd(E)=0A MeV, however, cannot be used at 



very low energies, since the DR cross sections are expected to be extremely small in that 
energy region. The systematics of the data suggest that a D may become essentially zero 
for E <10 MeV. We thus assume that Wd(E) increases linearly from zero at 10 MeV 
to the value of 0.4 MeV at £7=14.0 MeV. The strength W D {E) and the radius r D {E) 
parameters in the entire energy range E can then be rewritten as 



W D (E) 



0.0 

0.1(^-10.0) 

0.40 

0.40 



, r D (E) 



1.730 
1.730 

1.730- 0.03(^-14.0) 
1.508 



E < 10.0, 
10.0 < E < 14.0, 
14.0 < E < 21.4, 
21.4 < E. 



Together with a D =1.25 fm, Eq. (11) now defines W D (r; E) in the whole range of E. 

In generating the real part of the DR potential, Vd(t;E), by using the dispersion 
relation, we introduce an additional simplification of approximating the energy depen- 
dence of Wd{t;E) between £=14.0 and 21.4 MeV, where r D (E) changes with E. We 
assume a quadratic function of E for Wo(r;E) at each radial point r; Wr>(r;E) = 
a + b(E — E b ) + c(E — E b ) 2 , where a, 6, and c depend on r. We have confirmed that the 
approximation is accurate. Once this is done, the integration over E involved in Eq. (4) 
can be carried out analytically and one can get a closed form of Vo(r; E), 



V D (r;E) = V s (r; E s ) + -W D (r; E b )[e b ln\e b \ - e a ln\e a \] 

71 

+ -(W D (r; E c ) - W D (r; E b ))[e'Jn\<\ ~ 

7T 

-(W D (r; E c ) + W D (r; E b ) - 2W D (r; E m ))\e' c e' b (ln\e' c \ - ln\e' b \) + e' 6 ], 

71 

where e a and e b are the same as defined in Eq. (10) and 



(E - E b ) 



, and e' c = 



(E - E c ) 



(E c — E b ) ' (E c — E b ) ' 

with £ a =10.0 MeV, £ 6 =14.0 MeV , E C =2\A MeV, and E m = (E b + E c )/2. 



(12) 



(13) 
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Using the polarization potentials thus generated we perform the final calculations for 
elastic scattering, total DR and fusion cross sections and present the results in Figs. 4 and 
5, in comparison with the data. The data are fairly well reproduced by the calculations. 
The final calculated cross sections are essentially the same as those obtained in the initial 
X 2 -analysis, showing that the calculated cross sections do not sensitively depend on the 
real polarization potential, as we assumed in carrying out the x 2 analysis. We note that 
the fits to the elastic scattering and reaction cross sections (sum of the DR and fusion 
cross sections) are essentially the same as those obtained in Ref. ||. The fit to the 
elastic scattering data at the lowest energy £'=14.3 MeV is the worst among those shown 
in Fig. 4, but can be improved if we carry out a x 2 - ana lysis including only the elastic 
scattering data as the data to be reproduced. We made such an analysis, finding that the 
data were very well reproduced with r£>=1.93 fm, much larger than r^=1.72 fm obtained 
earlier. The DR cross section calculated with this r^i=1.93 fm, however, turned out to 
be <td=540 mb, about 3 times larger than the experimental value. This implies that one 
cannot improve the simultaneous fit to both the elastic and DR data any further. 

In summary, we have carried out simultaneous x 2 - ana lyses of elastic scattering, DR 
(breakup), and fusion cross sections for the 6 He+ 209 Bi system at near- Coulomb-barrier 
energies within the framework of an optical model that introduces two types of imaginary 
potentials, for DR and fusion, and determined the parameters of these potentials. The 
results indicate that the fusion potential exhibits a threshold anomaly very similar to 
that observed for tightly bound projectiles, but the magnitude at around the strong 
absorption radius is much smaller than the imaginary DR potential that does not show 
such an anomaly. Therefore, the resulting total imaginary potential does not show the 
anomaly. 
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FIGURE CAPTIONS 



Fig. 1. The experimental elastic and DR probabilities, P E and P D , respectively, as 
a function of the reduced distance d for the 6 He+ 209 Bi system at E cm = 18.5 and 21.9 
MeV. The data are taken from Ref. 9. The thin lines connecting Pi (i = E and D) values 
are only to guide the eye. 

Fig. 2. The values of Wf(E) (upper panel) and ro{E) (lower panel) extracted from 
the x 2 -analyses are plotted by the open and the solid circles, respectively. The solid lines 
represent Eqs. (7) and (8). 

Fig. 3. The real (upper panel) and the imaginary (lower panel) parts of fusion (dotted 
line) and DR (dashed line) potentials as functions of E at the strong absorption radius 
r = R sa =13.0 fm. The sum of fusion and DR potentials is plotted by the solid lines. 
The real parts of the potentials are calculated from Eq. (9) for Vf{E), and Eq. (12) for 
V D (E). 

Fig. 4. The ratios of the elastic scattering cross sections to Rutherford cross sections, 
calculated with our final optical potential for the 6 He+ 209 Bi system in comparison with 
the experimental data. The data are taken from Ref. 3. 



12 



Fig. 5. The calculated direct reaction and fusion cross sections with our final optical 
potential for the 6 He+ 209 Bi system in comparison with the experimental data. The data 
are taken from Refs. 3 and 15. 
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